Superfluidity of p-wacve and s-wave atomic Fermi gases in optical lattices 
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We consider p-wave pairing of single hyperfine state and s-wave pairing of two hyperfine states 
ultracold atomic gases trapped in quasi-two-dimensional optical lattices. First, we analyse superfluid 
properties of p-wave and s-wave symmetries in the strictly weak coupling BCS regime where we 
discuss the order parameter, chemical potential, critical temperature, atomic compressibility and 
superfluid density as a function of filling factor for tetragonal and orthorhombic optical lattices. 
Second, we analyse superfluid properties of p-wave and s-wave superfluids in the evolution from 
BCS to BEC regime at low temperatures (T ~ 0), where we discuss the changes in the quasi-particle 
excitation spectrum, chemical potential, atomic compressibility, Cooper pair size and momentum 
distribution as a function of filling factor and interaction strength for tetragonal and orthorhombic 
optical lattices. 
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I. INTRODUCTION 



Tunable optical lattices have been extensively used to 
study phase transitions in bosonic atomic gasesi^, since 
they allow the controlled manipulation of the particle 
density n, and of the ratio between the particle trans- 
fer energy t, and the interparticle interaction strength 
y i 3 j 4_ This kind of control is not fully present in standard 
fermionic condensed matter systems, and has hindered 
the development of experiments that could probe sys- 
tematically the effects of strong correlations as a function 
of n and t/V . However, fermionic atomic gases like 6 Li 
and K have been succesfully trapped, and their normal 
state and superfluid properties are beginning to be stud- 
iec&Si2iS. Because of the greater tunability of experi- 
mental parameters, novel superfluid phases may be more 
easily accessible in the experiments involving ultracold 
atomic gases. For instance, single hyperfine state (SHS) 
ultracold atomic systems are ideal candidates for the ob- 
servation of novel triplet superfluid phases and for testing 
theoretical models that were proposed earlier. Thus, it 
is only natural to propose that optical lattices could be 
used to study the normal state and superfluid properties 
of ultracold fermionic systems as a function of n, t/V 
and lattice symmetry. These systems are of broad inter- 
est not only for the atomic physics community but also 
for the nuclear, condensed matter and more generally for 
the many-body physics communities, where models for 
superfluidity have been investigated in various contexts. 

The interaction between induced dipole moments of 
atoms and the electric field of laser beams is used to trap 
atoms, particularly alkali atoms, in optical lattices. Al- 
kali atoms have only one electron (S — 1/2) out of closed 
shells. This electron is in a zero orbital angular momen- 
tum (L = 0) channel, and its total angular momentum 
J = L + S gives J = 1/2. The nuclear angular mo- 
mentum I and electron angular momentum J are com- 
bined in a hyperfine state with total angular momentum 
F = I + J which gives F = I ± 1/2 for alkalis. Fur- 
thermore, the electron and nuclear spins are coupled by 



the hyperfine interaction that splits the atomic levels in 
the absense of magnetic field Hm oc I • J. A weak mag- 
netic field causes Zeeman splitting of the hyperfine levels 
\F,rriF > with different mp, which can be made to corre- 
spond to pseudo-spin labels. Trapping neii*i2*i&±iii£ii& J 
twoiLiSiASiSSi, threeSi and four— hyperfine states were 
considered by several authors. However, it is also known 
that trapping more hyperfine states increases the num- 
ber of channels through which the gas can decay. There- 
fore, trapping one or two hyperfine states of ultracold 
fermionic gases is experimentally more plausible. 

For paired identical fermions, the Pauli exclusion prin- 
ciple requires the total pair wave function to be anti- 
symmetric. The total orbital angular momentum should 
be odd for pseudo-spin symmetric pairs and even for 
pseudo-spin anti-symmetric ones. Therefore in the case 
of trapping two hyperfine states (THS), s-wave scatter- 
ing of atoms between fermions from different hyperfine 
states is dominant. Thus, one expects that the super- 
fluid ground state of such two-component Fermi gases 
to be s-wave and pseudo-spin singlet. Presently there 
is experimental evidence that 40 K (Refi2i24) and 6 Li 
(Refi 25 i 26 i 27 i 28 i 29 ) can form weakly and tightly bound 
atom pairs, when the magnetic field is swept through 
an s-wave Feshbach resonance. 

However, the properties of SHS ultracold fermions 
and their possible superfluid behaviour are beginning to 
be investigated22i2ii22i2&2i. These systems are probably 
the next frontier for experiments with ultracold atoms. 
When identical fermionic atoms are trapped in a single 
hyperfine state, the interaction between them is strongly 
influenced by the Pauli exclusion principle, which pro- 
hibits s-wave scattering of atoms in identical pseudo-spin 
states. As a result, in SHS degenerate Fermi gases, two 
fermions can interact with each other at best via p-wave 
scattering. Thus, one expects that the superfluid ground 
state of such SHS Fermi gases to be p-wave and pseudo- 
spin triplet. 

In the p-wave channel, if the atom-atom interactions 
are effectively attractive then the onset for the forma- 
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tion of Cooper-pairs in three dimensions occurs at a 
temperature!^ T c « e F exp[-7r/2(fc F a sc ) 3 ] in the BCS 
regime, where €f is the Fermi energy and a sc is the p- 
wave scattering length. Unfortunately, this temperature 
is too low to be observed experimentally. However, in 
the presence of Feshbach resonancesifi^, p-wave interac- 
tions can be enhanced, and the critical temperature for 
superfluidity is expected to increase to experimentally 
accessible values. On the other hand, we show that the 
pseudo-spin triplet (p-wave) weak coupling limit in op- 
tical lattices (like in the singlet case£) may be sufficient 
to produce a superfluid critical temperature that is ac- 
cessible experimentally. Furthermore, several interesting 
superfluid properties can be investigated in p-wave sys- 
tems as the system is tuned from the BCS to the BEC 
rcgimea 10 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 ^ 7 ! 18 . In this paper, we address 
superfluid properties of ultracold fermionic atoms in op- 
tical lattices for both SHS p-wave and THS s-wave states 
as a function of atom filling factor, interaction strength, 
temperature and lattice symmetry. 

The rest of the paper is organized as follows. In SeclTTl 
we discuss the effective action method and the saddle 
point approximation for a lattice Hamiltonian with at- 
tractive interactions in the p-wave and s-wave channels. 
We also derive the order parameter, critical temperature 
and number equations in the saddle point approximation 
and discuss pseudo-spin triplet states in the d- vector for- 
malism in this section. In Sec. IIIII we analyse superfluid 
properties of SHS p-wave and THS s-wave symmetries in 
the strictly weak coupling BCS regime. There, we discuss 
the order parameter, chemical potential, critical temper- 
ature, atomic compressibility, and superfluid density as 
a function of filling factor for tetragonal and orthorhom- 
bic optical lattices. In Sec. IIV1 we analyse superfluid 
properties of SHS p-wave and s-wave superfluids in the 
evolution from BCS to BEC regime at low temperatures 
(T w 0). There, we discuss changes in quasi-particle ex- 
citation spectrum, chemical potential, atomic compress- 
ibility, Cooper pair size and momentum distribution as 
a function of filling factor and interaction strength for 
tetragonal and orthorhombic optical lattices. Finally a 
summary of our conclusions is given in Sec. [V] 



II. LATTICE HAMILTONIAN 

In this manuscript, we consider quasi- two-dimensional 
optical lattices with a periodic trapping potential of 
the form U(r) = i Uo,i cos 2 (kiXi), with Uo, z ^> 
min{J7o,a:, C^o,j/}i which strongly suppresses tunneling 
along the z direction. This is a non-essential assumption, 
which just simplifies the calculations, but still describes 
an experimentally relevant situation. Here X{ — x, y, or z 
labels the spatial coordinates, fcj = 2ir/Xi is the wave- 
length, and Uo,i is the potential well depth along direction 
Xj, respectively. The parameters C/o,i are proportional to 
the laser intensity along each direction, and it is typically 
several times the one photon recoil energy Er such that 



tunneling is small and the tight-binding approximation 
can be used. 

Thus, in the presence of magnetic field h, we consider 
the following quasi-two-dimensional lattice Hamiltonian 
(already in momentum space) for an SHS p-wave Fermi 
gas 

H = J2 C(k)4 T «kT + \ £ V p (k, k')< q & k ,, q , (1) 

k k,k',q 

where the pseudo-spin | labels the trapped hyper- 
fine state represented by the creation operator ojL-, 

and b l.q = a k+q/2,T a -k+q/2,T- FurtnCmlOTC > £0) = 

e(k) — jj, describes the tight-binding dispersion e(k) = 
—t x cos(k x a x ) — t y cos(k y a v ) — t z cos(k z a z ), with jl = 
H + g^Bh — VH- Here, U is a tranfer energy along the i-th 
direction, /i is the chemical potential, gps/i is a Zeeman 
energy, Vh is a possible Hartree energy shift. We also 
assume that mhi{t x ,ty} 3> t z , and that V(k, k') is the 
pseudo-spin triplet pairing interaction between fermions. 

In the nearest neighbour approximation, the lattice in- 
teraction in the single hyperfine state case has only a 
p-wave (triplet) component, which is given by 

V p {k, k') = -2 V o,i sin{kai) sin(^a 4 ), (2) 

i— x,y 

where Vo,i > is the effective interaction strength and 
at is the corresponding lattice length along the i th direc- 
tion. Notice that Eq. J3J) has the necessary symmetry 
under the Parity operation, where either k — ► — k or sim- 
ilarly k' — >■ — k' leads to — V^(k, k'). Furthermore, the 
interaction V^(k, k') is invariant under the transforma- 
tion (k, k') — > (— k, — k'). This is a necessary condition 
for an SHS p-wave interaction since the interaction char- 
acterizes a triplet channel, and it has to reflect the Pauli 
principle. This interaction can be written in a separable 
form as 

7 p (k,k') = rt(k)vr(k') (3) 

where the interaction strength matrix has elements 
Vij = —2Vo t iSij and the symmetry vector is T^(k) = 
[sin(k x a x ),sin(kya y )]. Notice that, the interaction 
Vp(k, k') is not separable as a product of scalar func- 
tions of k and k' in the usual sense, but it is separable 
in terms of vector functions. 

In addition, we compare SHS p-wave superfluids with 
THS s-wave superfluids. In the nearest neighbour ap- 
proximation, the lattice interaction term in the THS case 
leads to singlet pairing channels s, extended-s and d-wave 
terms, and a separate triplet pairing THS p-wave chan- 
nel. However, the s-wave channel is expected to dominate 
in the absence of any exotic mechanism as it corresponds 
to the case of minimal free energy. The Hamiltonian for 
a THS s-wave Fermi gas is given by 

F = ]T£(k)4 CT a kCT + i ]T Kbt^,^ (4) 

k.tr k,k' ,cr,q 
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where the pseudo-spins a — (T, 4) label the trapped hy- 
perfine states represented by the creation operator a ka , 

aild & Lq = «lc+q/2^ at -k+q/2,-<T- Here > V » = ~\ V oA ^ 

the attractive s-wave interaction. 

In the following section, we describe the effective ac- 
tion method for the SHS p-wave Hamiltonian in Eq. 
and we discuss its saddle point order parameter and num- 
ber equations, as well as the critical temperature. This 
method can be applied to the s-wave Hamiltonian in 
Eq. |0J, and similar equations can be derived for s-wave 
case, but we do not repeat this derivation here. 



Effective Action Formalism 



In the imaginary-time functional integration formal- 



lsrnS 

function can be written as 



(/3 = l/T and units % = fee = 1) ; the partition 



Z = / D\a\a]e 



with an action given by 



(5) 



(6) 



The Hamiltonian given in Eq. Q can be rewritten in the 
form 

H(r) =^^(k)a| cT (r)a kT (r)+^Bt(r)^B q (T) (7) 



has r-independent and T-dependent parts where Aj = 
( A o,*, A o, y ) and At (<?) = i A x(q), A y (q% respectively. 

Performing an expansion in S c s to quadratic order in 
the vector A, we obtain 



5', 



— Sn 



i^At^F- 1 ^), 



(11) 



where So is the saddle point action, the 4-component vec- 
tor At(g) is such that A!(q) = [A^(q),A(-q)} and P _1 (g) 
is the inverse fluctuation propagator. The fluctuation 
term in the action leads to a correction to the thermo- 
dynamic potential, which can be written as fi gaU ss = 
n + fi fluct with fifluct = f3- 1 Y, q ^det[F- 1 (q)/2}. In 
weak coupling for all temperatures and for all couplings 
at low temperatures, it can be shown that^S Gaussian 
corrections Sfi uct to the gaussian action S gauss are small 
in the quasi-two-dimensional case, and therefore, we con- 
sider only the saddle point action 



S = 2pAlv- 1 A + Y / 



£(k) 



Trln-Go" 

2 



,(12) 



where the inverse Nambu propagator becomes 

G^ 1 = iw e a - £(k)a 3 + A o r(k) ( x_ + rt(k)A a+. (13) 

Notice that, this will not be the case close to the critical 
temperatures T c for intermediate and strong interactions, 
since the inclusion of fluctuation corrections change the 
number equation considerably and are necessary to pro- 
duce correct behaviour—. 



with B q (r) = J2k r(k)&k. q ( T )- We then introduce 
Nambu spinor V^(p) = { a v -\i a -pt)> where we use p = 
(k, wg ) to denote both momentum and fermionic Mat- 
subara frequency wi = (21 + l)7r//3 and use a Hubbard- 
Stratonovich transformation to decouple fermionic and 
bosonic degrees of freedoms. Performing the integration 
over the fermionic part (D[ip> , tp]) leads to 



5, 



c rr 



E 



^(k)^-Trln^G- 1 



(8) 



where we use q = (q, vg ) with bosonic Matsubara fre- 
quency vg = 2£w//3 and define the bosonic vector field 
$>(p-p'). Here, 



G 1 = *\ q) T { ^f)a_ + T\ P -±fm- q )a, 
+ [iw e cr - £(k)<T 3 ] 6 p y 



(9) 



is the inverse Nambu propagator and a± — (o\ ± er 2 ) /2 
and Ci is the Pauli spin matrix. The bosonic vector field 



= Ao<5 g ,o + A(g) 



(10) 



B. Saddle- Point Equations 

The saddle point condition 5Sq/SA^ = leads to a 
matrix equation for the order parameter 

A = MA , (14) 

where matrix M has the following matrix elements 

Mij = J2 VQ/l g^^pte) tanh ggW . (15) 

Here, we introduce the quasi-particle energy 

£(k) = (£ 2 (k) + |A(k)| 2 )^ (16) 

and the scalar order parameter 

A(k)=rt(k)A 0) (17) 

which is naturally separable in temperature T and 
momentum k. The critical temperature, T c — 
max{T c a; , T c y \, is determined from the condition 
dot M = I in Eq. I|14|) . and can be written as 



n u-^E 



Sin2(fc4fll) tanh^ k) 



2T C . 



(18) 
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Both the order parameter and critical temperature 
equations have to be solved simultaneously with the num- 
ber equation N — —d£l/dfi where is the full thermo- 
dynamic potential. This leads to two contributions to 
the number equation N w N gauss = Nq + Nft uct . Here, 
Nq = —dflo/d^i is the saddle point number equation, 
where flo — So/f3 is the saddle point thermodynamic 
potential, and iVfl uc t = — 9f2fl uct /9/x is the fluctuation 
contribution. In the low temperature limit (T rs 0) for 
any coupling, or in the weak coupling limit for any tem- 
perature (T < T c ) the fluctuation contribution 



fl[detF- 1 (g)] 



1 



detP^fg) 



(19) 



is small and negligible compared to the saddle point value 
iV(j2£. Thus, the number equation becomes N w N Q , 
given by 



N 



(20) 



where n(k) is the momentum distribution defined by 



ra(k) 



E(k) 2 



(21) 



For comparison, we discuss next the s-wave case. Us- 
ing the s-wave Hamiltonian in Eq. J3J, we follow a sim- 
ilar procedure and derive the order parameter, critical 
temperature and number equations in the saddle point 
approximation. The order parameter equation is given 
by 



A Q , S = V 0}S ^ 



A 



0,s 



2E{k) 



tanh 



(22) 



where E(k) — y/£ 2 (k) + |Ao jS | 2 is the quasi-particle en- 
ergy spectrum. The critical temperature is determined 
when Aq jS = 0, which leads to 



l = £^tanh^ 



2£(k) 



2T r 



(23) 



Both order parameter and critical temperature equations 
have to be solved simultaneously with the number equa- 
tion 



N 



n(k), 



(24) 



where the momentum distribution n(k) is defined in 
Eq. HU). 



C. Spin Triplet (SHS p-wave) States 



In general, the pseudo-spin triplet order parameter can 
be written in the standard form 36 



O(k) 



-di(k)+tda(k) d 3 (k) 
d 3 (k) d 1 (k) + id 2 (k) 



(25) 



where d(k) is an odd vector function of k. In our SHS 
p-wave interaction Off (k) = A(k), therefore, d 3 (k) = 
and di(k) = — id2(k) which leads to d(k) = di(k)(l, i, 0). 
Within the irreducible representations of the T>4h (D2/1) 
group in the tetragonal (orthorhombic) lattices, — our 
exotic SHS p-wave state corresponds to the 3 E u (n) rep- 
resentation with a d- vector given by 



d(k) = /(k)(l,i,0), 



(26) 



where /(k) = AX + BY, and X and Y are sm(k x a x ) 
and sin(k y a y ), respectively. Notice that, this state also 
breaks time-reversal symmetry, as expected from a fully 
spin-polarized state. 

In the tetragonal lattice, the stable solution for our 
model corresponds to the case A = B 7^ 0, and thus 
to the 3 E u (d) representation, where spin-orbit symme- 
try is preserved, but both spin and orbit symmetries are 
independently broken. In the orthorhombic lattice, the 
stable solutions correspond to either A 7^ 0, B = or 
A = 0, B 7^ 0, thus leading to the 3 E u (b) representation. 
Notice that, depending on the lattice anisotropy, there 
are three distinct solutions to the order parameter equa- 
tion Eq. i|14|) in relation to Eq. I|26|) . The first solution 
(case 1) occurs when Aq^ x 7^ and Ao,j, = 0, which cor- 
responds to a d-vector with A 7^ and B = 0. Similarly, 
a second solution (case 2) occurs when A 0jX = and 
Ao.j, 7^ 0, and it corresponds to a d-vector with A = 
and B ^ 0). Third and final solution (case 3) occurs 
when Aq iX = Ao, a 7^ 0, which corresponds to a d-vector 
with A = B 7^ 0. In a tetragonal lattice both directions 
are degenerate (case 3), but even small anisotropics in 
the optical lattice spacings, transfer energies (t x ,t y ), or 
optical lattice potentials (Uo yX , Uo tV ) lifts the degeneracy 
and throws the system into either case 1 or 2. 



III. SUPERFLUID PROPERTIES IN THE 
WEAK COUPLING BCS REGIME 



Our main interest is in tetragonal (square) lattices, 
however, we also want to investigate the effects of small 
anisotropies in optical lattice lengths. For this purpose 
we investigate five different cases that may be encoun- 
tered experimentally 

Case (I) corresponds to lattice spacings a x = a y , to 
transfer integrals t x = t y , and to interaction strengths 
Vq iX = Vo tV . Case (II) corresponds to a x = a y , t x 7^ t y 
and Vo, x — Vo tV . Case (III) corresponds to a x = a y , 
t x 7^ t y and V~o, x 7^ Vq iV . Case (IV) corresponds to a x 7^ 
a y , t x 7^ t y and Vq^ x = Vo, y . Case (V) corresponds to 
a x 7^ a y , t x 7^ t y and Vq. x 7^ Vo. y . Notice that cases 
(II) and (IV) are equivalent except for a unit cell volume 
normalization. The same is true for cases (III) and (V). 
Case (I) is effectively a tetragonal lattice, while cases (II) 
through (V) are effectively orthorhombic lattices. Notice 
that even though the optical lattice spacings a x — a y in 
cases (II) and (III) , these cases correspond effectively to 
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" orthorhombic" lattices because the transfer energies t x 
and t y are different. 

Experimentally it may be easier to use only one type 
of laser with specified wavelength, and thus generate a 
lattice with a x — a y , and change the focus width in one 
direction (say along x) while keeping the width in the 
other direction fixed (say along y). The change in focus 
width along the x direction will modify the transfer inte- 
gral t x by either reducing it or enlarging it with respect 
to t y . This type of experiment can cover cases (I), (II) 
and (III). A second type of experiment (a bit harder) 
could use two types of lasers with different wavelengths 
and produce a lattice with a x ^ a y . Further control of 
the focal widths can produce the situations encountered 
in cases (IV) and (V). 

We concentrate here on cases (I), (II) and (III), where 
we set a x — a y = a. For case (I) we choose t x = t y = t, 
and Vo,x = Vo,y — Vq. For case (II) we choose t x = 
(1 + a)t, t y = t, and Vo, x — Vo, v — Vq. For case (III) 
we choose t x = (1 + a)t, t y = t, and V$, x — (1 + ct)Vo, 
Vo.y — Vq. Here we take Vq = 0.6_E with E = 2t. 
In these cases, the transfer energies and the parame- 
ters a can be determined by using exponentially decay- 
ing on-site Wannier functions and the WKB approxima- 
tion. When the lattice spacings a x = a y = a are the 
same but the optical lattice potentials Uo, x > Uq jV = Uq 
are slightly different, the transfer integral t x /t is propor- 
tional to exp [— tv(Uq.x — Uq)/ WUo~xEr], in the limit of 
Uq,x — Uq <C Uq, x - Here, E R is the one photon recoil 
energy, while Uq and Uo tX are defined with respect to E g , 
which is the energy of the lowest reference state of the 
trapping potential. In the following sections, we discuss 
two possibilities where Uo, x — Uq = 0.05-Er correspond- 
ing to a « 0.1 and Uq x — Uq = 0.2E R corresponding to 
a « 0.6 for U « 2.5E R . 

For cases (IV) and (V) we set a x = a(l — 8) and 
a y = a, with < S <C 1, respectively. For case (IV) 
we choose t x = (1 + a)t, t y = t, and Vq, x = Vq iV = Vq. 
For case (V) we choose t x = (1 + a)t, t y — t, and 
Vq. x = (l+a)Vo, Vo tV — Vq. In these cases, the parameter 
a can also be determined by using exponentially decaying 
on-site Wannier functions and the WKB approximation. 
When the optical lattice potentials Uq <x = Uo tV = Uq 
are the same but the lattice spacings a x < a y = a 
are different, the transfer integral t x /t is proportional 
to exp (—2ttS^Uq/E r ), for small 6. In the following sec- 
tions, we discuss two possibilities where S = 0.01 cor- 
responding to a w 0.1 and 6 = 0.05 corresponding to 
a w 0.6 for Uq fts 2.5E R . As we increase (decrease) the 
ratio <5, both interactions and tunneling rate increase (de- 
crease) in x direction. 

Notice that a different choice of on-site Wannier func- 
tions does not change our general conclusions, the only 
qualitative difference is that t x and Vq iX are different 
functions of a. 



A. Density of Fermionic States 

The critical temperature and superfluid properties of 
Fermi systems depend highly on the order parameter 
symmetry, as can be seen by rewriting Eq. (|18f) as 

„ r t*+t tanh|^ 
1 = II V 0li / de J Am(e), (27) 

i—x,y —tz—t \ r 6 / 

where we define an effective density of states (EDOS) 

d pAs) = ^<5[e-£(k)][V2sin(M J )] 2 . (28) 

k 

Here y/2sm(kia,i) is the symmetry factor related to the 
SHS p-wave order parameter. We transformed discrete 
summations over k-space to continuous integrations to 
obtain 

D VtX (e)= f ^[tl-{e-tco S {k v a)f] 1J \ (29) 

Ji I* l x 

where the integration region is restricted by 

| £ -tcos(fc y a) | ^ 1 ^ ^ firgt Brillouin zone (1BZ). 

Plots of EDOS are shown in Figs. ^ an( i f° r three 
cases: a = (hollow squares), a — 0.1 (solid squares), 
and a — 0.6 (line). 

In Fig. ^i, we plot EDOS which is smooth and maxi- 
mal at half-filling in square lattices (a = 0) . For the SHS 
p-wave symmetry discussed, Vo t iD Pi i(e) plays the role of 
a dimcnsionlcss coupling parameter which controls the 
critical temperature. This is simply because it is much 
easier to form Cooper pairs with a small attractive inter- 
action (and lower the free energy of the Fermi system) 
when the density of single fermion states is high. Ad- 
ditionally, EDOS decreases and finally vanishes at the 
band edges where a small ratio of T c /Tp was predicted 
from continuum models. However, this is not the case 
around half-filling, and we expect weak interactions to 
be sufficient for the observation of superfluidity. 

In contrast, the symmetry factor is 1 in the s-wave case 
and EDOS becomes 

D( £ ) = ^{[ E -e(k]], (30) 

k 

which is identical to the density of single fermion states 
(DOS) of normal fermions. A similar calculation yields 

D(s) [tl ~ (e ~ t cos(fc,a)) 2 ] ~ 1/2 , (31) 

where the integration region is again restricted by 

| £ -tcos(fc y a) | ^ 1 ^ ^ plotg of ^ ghown 

in Fig. ^ an d Fig. ^1 for three cases: a — (hollow 
square), a = 0.1 (solid square), and a = 0.6 (line). 

For instance, notice that in the tetragonal case of a two 
dimensional lattice, the DOS at half-filling is very large 
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(b) 



(a) 



(b) 



0.5 



0.016 





FIG. 1: Plots of ED OS versus reduced energy e r = e/Eo for 
SHS p-wave symmetry when a) a = (hollow square), and b) 
a = 0.1 (solid square) and a — 0.6 (line). In addition plots 
of DOS for s-wave symmetry when c) a = (hollow square), 
and d) a = 0.1 (solid square) and a — 0.6 (line) are shown. 



and tends to infinity logarithmically (Fig. ^;). However, 
away from half-filling the DOS decreases rapidly to a con- 
stant at the minimum and maximum band edges. This 
constant DOS for very low or very high filling factors is 
equivalent to the DOS of the continuous (homogeneous) 
systems since at the bottom or at the top of the band, the 
single fermion dispersion energy can be approximated by 
a parabola in momentum space. Furthermore, for s-wave 
superfluids, VqD(e) is the dimensionless coupling param- 
eter which controls the value of critical temperature T c . 
This important quantity (T c ) is discussed next. 



B. Critical Temperature 

Plots of reduced critical temperature T r = T c /Eq as a 
function of the number of atoms per unit cell (N c ) are 
shown in Fig. [21 for SHS p-wave and s-wave superfluids 
in cases (I), (II), and (III). Here, T c is the critical tem- 
perature calculated from Eqs. I|18fl or (|27|l and Eq = 2t 
is the half-filling Fermi energy for the square lattice with 
respect to the bottom of the band. 

The square lattice case (I) is shown in Figs. [2K and [2: 
(hollow squares) for SHS p-wave and s-wave superfluids, 
respectively. Notice that T r is maximal at half-filling, 
and that it has a value 0.01, which is much higher than 
the theoretically predicted T c from a continuum model, 
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FIG. 2: Plots of reduced critical temperature T r = T c /Eo 
versus filling factor N c for SHS p-wave symmetry are shown 
for cases a) (I) and (II) and b) (I) and (III). In addition, plots 
for s-wave symmetry is shown for cases c) (I) and (II), and 
d) (I) and (III). Plots for a — 0,0.1 and 0.6 are shown with 
hollow squares, solid squares and line, respectively. 



and comparable to experimentally attainable tempera- 
tures T/T F « 0.042&. This implies that the superfluid 
regime of SHS fermion gases may be observed experimen- 
tally in a lattice, even in the limit of weak interactions 
(BCS regime). The observability of a superfluid tran- 
sition in SHS p-wave Fermi systems is clearly enhanced 
when the system is driven through a Feshbach resonance 
(in a lattice or in the continuum), as T c is expected to 
increase further in this case, however our calculations 
indicate that the weak interaction (BCS) limit may be 
sufficient in the lattice case. This possibility for Fermi 
gases in optical lattices has a parallel in the experimental 
observation of superfluid-insulator transitions in weakly 
interacting Bose systems in optical lattices^. 

In case (II), while we keep interaction strengths Vq iX — 
Vo,y = Vq the same, we vary only the transfer energy 
along x direction t x = (1 + a)t. This case is shown in 
Fig. [21i and Fig. [2J; for SHS p-wave and s-wave superflu- 
ids, respectively. Here, solid squares (line) correspond to 
a = 0.1 (a = 0.6). 

For SHS p-wave superfluids, an infinitesimally small 
anisotropy in the transfer energy ratio t x /t y causes a dis- 
continous jump in EDOS at half-filling because of the 
symmetry factor enhancement, however T c at half-filling 
is the same as in case (I), because the discontinuity at 
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EDOS is limited to a single point. As the anisotropy ra- 
tio increases (a = 0.1), the discontinuous point in EDOS 
at half-filling expands into a region where EDOS is higher 
than in case (I). Thus over a narrow range around half- 
filling, T c is larger than in case (I), but varies smoothly 
as a function of a or N c . Further increase of a reduces 
EDOS near half-filling pushing T c down. This behavior is 
characteristic of the SHS p-wave state through its symme- 
try factors. In contrast, for s-wave superfluids shown in 
Fig-Gt, T c decreases for fixed N c around half-filling with 
increasing anisotropy (at least for < a < 0.6). Fur- 
thermore, T c develops two maxima which are symmetric 
around half-filling. This behavior in T c is also related 
to the DOS of the normal fermions which is shown in 
Fig. ^ Notice here that, when there is an anisotropy in 
the transfer energy, T c near half-filling corresponding to 
an SHS p-wave pairing is considerably higher than the T c 
of an s-wave pairing for the same parameters of interest 
(assuming that Vo, s = Vo tP = Vq). 

In case (III), we vary both the interaction strength 
Vq, x = (1 + a)Vo and the transfer energy t x = (1 + a)t 
along £ direction. This case is shown in Figs. [2b and[2Ji 
for SHS p-wave and s-wave superfluids, respectively. In 
these figures it is assumed that Vo, s = Vq iP = Vb, x - In 
case (III) , both SHS p-wave and s-wave superfluids have 
critical temperatures that are much larger than in case 
(II), because of the increase in Vq, x , in addition to the 
density of states effect, discussed above. Notice again 
that T c near half-filling corresponding to an SHS p-wave 
pairing is also considerably higher than T c for an s-wave 
pairing for the same interaction parameters. 

In this section, we established that the critical temper- 
atures for SHS p-wave superfluids in optical lattices can 
be comparable or larger to the critical temperatures of 
s-wave superfluids, depending on the lattice anisotropy 
and interaction strength. Thus, we expect that SHS p- 
wave superfluids may be experimentally attainable. In 
the next sections, we discuss some experimentally rel- 
evant observables which could be used to identify the 
superfluid phase of our exotic SHS p-wave state. 



C. Order Parameter and Chemical Potential 

In this section, we discuss the low-temperature be- 
havior of the order parameter A(k) = A x sin(k x a x ) + 



(a) 



Ao j?/ sin(fcj,a a ) defined in Eq. 
state. In case (I) A Q ^ X = A 



(jH|) for the SHS p-w&ve 
9 ^ 0, while in cases (II) 



and (III) Ao,i 7^ and Ao, a = 0. Here, we also discuss for 
comparison the singlet s-wave case where A(k) = Ao jS , 
with A , s ^ for cases (I), (II) and (III). The re- 
duced amplitudes of the SHS p-wave order parameter 
A,. = Ao, x /Eq versus N c are shown in Figs.|3^, and^- In 
addition, we plot s-wave A,. = Aq. s /Eq in Figs.[2t and[3Ji 
for cases (I), (II), and (III). Notice that, the qualitative 
behaviour of A r as a function of N c is very different from 
that of the chemical potential [i r = fti/ Eq, which is dis- 
cussed next. 
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FIG. 3: Plots of zero temperature reduced order parameter 
amplitude A r = Aq <x /Eo versus filling factor iV c for SHS 
p-wave symmetry for cases a) (I) and (II) and b) (I) and 
(III). In addition, plots for s-wave symmetry A r = Aq iS /Eo 
is shown for cases c) (I) and (II), and d) (I) and (III). Plots 
for a — 0,0.1 and 0.6 are shown with hollow squares, solid 
squares and line, respectively. 



Plots of the low-temperature /i r versus N c are shown 
in Fig. 2] for SHS p-wave and s-wave superfluids, respec- 
tively. Here, we present only the corresponding chem- 
ical potentials for cases (I) (hollow squares) and (II) 
(solid squares), where a = and a — 0.1, respec- 
tively. Chemical potentials for case (III) are very sim- 
ilar to those of case (II). Notice that ft, is always within 
the limits of the energy dispersion of the optical lattice, 
— t x — t < ft < t x + t, which characterizes the weak cou- 
pling (BCS) limit. 

Although the chemical potentials of SHS p-wave and 
the singlet s-wave cases look similar at first glance, they 
are not. The major qualitative differences between the 
chemical potentials of the SHS p-wave and the singlet 
s-wave cases are more clearly seen in the derivative 
dft/ dN c . The appropriate thermodynamic quantity that 
is related to this derivative is the atomic compressibility 
to be discussed next. 



D. Atomic Compressibility 

The isothermal atomic compressibility at finite tem- 
peratures is defined by kt{T) — —(dV/dP)T,N/V where 
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FIG. 4: Plots of zero temperature reduced chemical potential 
fi r = p./Eo versus filling factor N c for a) SHS p-wave, and b) 
s-wave superfluids. Plots for a — and 0.1 are shown with 
hollow and solid squares, respectively. 



V is the volume and P is the pressure of the gas. This 
can be rewritten as 



kt(T) 



1 



d 2 n 

dp 2 



T,V 



1 



giVc 



The partial derivative cW c / <9/i is 



giVc \ - 
dpi ^ 2E 3 (k) 



tanh « 



■Em 



r.v 



g 2 (k) 

£ 2 (k) 



(32) 



(33) 



where Y(k) = (/3/4)sech 2 [/3E(k)/2] is the Yoshida func- 
tion. 

Plots of the reduced isothermal atomic compressibility 
K r = kt(T)/kq are shown in Fig. for both SHS p-wave 
and s-wave superfluids for cases (I) and (II). For case 
(III), k t is both qualitatively and quantitatively similar 
to case (II), and it is not shown here. The normalization 
constant kq is the SHS p-wave isothermal compressibility 
evaluated at T = and half-filling. 

In case (I), K r has a peak at half- filling and low tem- 
peratures, and a hump at T c for SHS p-wave superflu- 
ids. Notice the strong temperature dependence of n r 
near half-filling. In the s-wave case there is only a hump 
both at T c and at low temperatures, and a very weak 
temperature dependence for all filling factors shown. In 
case (II), K r has the additional feature that the central 
peak (hump) splits into two due to degeneracy lifting, for 
the SHS p-wave case. Notice again the drastic difference 
between the values of K r for N c w 0.45 and N c ps 0.55 
for T and T « T c . However, in the s-wave case the 
original hump also splits into two, but there is very weak 
temperature dependence of K r . 

To understand the peaks of n r at T « 0, and its humps 
at T = T c limits for the SHS p-wave case, we indicate that 
dN c /dp [see Eq. (|33H ] has two contributions, where the 
first term dominates at T = 0, and only the second term 
survives at T = T c . Therefore, the first term is responsi- 
ble for the peaks, and the second term is responsible for 
the humps. 
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FIG. 5: Plots of reduced compressibility n r = kt(T)/ko 
versus filling factor N c for SHS p-wave symmetry in cases 
a) (I) and b) (II). In addition, plots for s-wave symmetry 
are shown in cases c) (I) and d) (II) at temperatures T — 
(squares), and T — T c (line). Plots for a — and 0.1 are 
shown with hollow and solid squares, respectively. 



The peaks at T w 0, can be understood by noting that 
the first contribution to kt(T) in Eq. 113311 can be written 
as 



k t {T = 0) 



1 ^ n(k)[l-n(k)] ^ (34) 



E(k) 



where n(k) is the momentum distribution and defined in 
Eq. (|21ll . Therefore, the peaks are due to non- vanishing 
n(k)[l — n(k)] in regions of k-space where E(\&) vanishes 
and n(k) is rapidly changing. 

In case (I), the integrand n(k)[l — n(k)]/£'(k) has 4 
k-space points (0,±7r), and (±7t,0) in the first Brillouin 
zone (1BZ) where it diverges only when the chemical po- 
tential p, — (N c — 0.5). Similarly in case (II), the in- 
tegrand diverges at 2 k-space points (±7r, 0) in the 1BZ 
when p = t x -t fts 0.021 (N c sa 0.55). Furthermore, 
the integrand diverges at 2 k-space points (0, ±7r) when 
p = -t x + t& -0.021 (N c w 0.45). For other values of p 
(other filling factors), the integrand is well-behaved and 
does not produce additional peaks in K r . 

In contrast, the compressibility peaks at T = do 
not exist for s-wave superfluids (see Figs. and [5JI), 
where the quasi-particle excitation spectrum -E(k) is al- 
ways gapped. Generally speaking, we expect compress- 
ibility peaks for nodal superfluids or superconductors, 



9 



where quasi-particle energy spectrum vanishes in regions 
of k-space. 

At T — T c , kt(T) is dominated by the second term of 
the integrand (see Eq. iJBEJ) 



/CT (T = T c ) = -J-^sech ; 



2T C N? 



(35) 



Therefore, the humps are not related to the order pa- 
rameter (since at T = T c , A(k) = 0), but are due to the 
peaks appearing in the single fermion (normal) DOS (see 
Fig. [Tt and Fig.QJi). This is simply because n T (T = T c ) 
can be written in terms of DOS (see Eq. 1301) for both 
SHS p-wave and s-wave symmetries. Notice that, while 
DOS have only one peak at half-filling in case (I) (lead- 
ing to one hump in K r ), DOS has two peaks in cases (II) 
and (III) (leading to two humps in n r ). Thus, at T = T c 
the humps correspond to a Van Hove phenomenon^, but 
the atomic compressibility is smooth. Furthermore, we 
expect humps at T — T c for all pairing symmetries, since 
these humps are related only to normal state properties 
but not to the symmetry of the pairing interaction. 

Theoretically, the calculation of isothermal atomic 
compressibility kt (T) is easier than the isentropic atomic 
compressibility Ks{T). However, performing measure- 
ments of ks{T) may be simpler in cold Fermi gases, since 
the gas expansion upon release from the trap is expected 
to be nearly isentropic. Fortunately, Ks(T) is related to 
kt(T) via the thermodynamic relation 



(36) 



where Kt(T) > ks(T) since specific heat capacitites 
Cp(T) > CV(T). Furthermore, at low temperatures 
(T fn 0) the ratio C P (T)/C V (T) « 7 is a constant, and 
therefore, ks(T w 0) cx kt{T ~ 0). Thus, we expect 
qualitatively similar behaviour in both the isentropic and 
isothermal compressibilities at low temperatures (T k 0), 
where appearence of such peaks could be used as a sig- 
nature of SHS p-wave superfluidity. 

The measurement of the atomic compressibility could 
also be performed via an analysis of particle density fluc- 
tuations. As it is well know from thermodynamics^ the 
isothermal atomic compressibility is connected to density 
fluctuations via the relation 



(nfT 
(V) 



«t(T), 



(37) 



where (n) and (V) are average density of atoms and vol- 
ume of the condensate, respectively. From the measure- 
ment of density fluctuations the isothermal compressibil- 
ity can be extracted at any temperature T. 

Furthermore, the spin susceptibility tensor component 
% W (T) can be written as a spin-spin response function 
to a magnetic field h = hfj applied along an arbitrary 
redirection^. In the case of SHS p-wave superfluids, 



d 2 n 2 2 dN c 



Therefore, for the SHS p-wave case, x w (T) is directly 
related to the isothermal atomic compressibility and is 
given by [N 2 /{gpB) 2 ]KT{T). Thus, we expect similar 
effects in the magnetic spin susceptibility as in the case 
of isothermal atomic compressibility discussed above. In 
cold atoms, the measurement of the spin susceptibility 
will be most likely achieved by using techniques similar 
to nuclear magnetic resonance (NMR), where Xvv(.T) can 
be measured via the Knight shifta£. 

Specifically, while we expect only one peak at exactly 
half- filling in the tetragonal lattices at T = 0, two peaks 
will appear in the orthorhombic lattices for a SHS p- 
wave superfluid. However at T = T c , these peaks dis- 
appear and turn into humps. Notice that the filling fac- 
tor dependence of X W (T) will be qualitatively different 
from kt(T). This is because while Xnni^) IS symmetric 
around half-filling (N c = 0.5), Kt(T) is not. The relation 
between the particle compressibility and the spin suscep- 
tibility given above is not valid for the s-wave case. Their 
relationship is more complicated in this case, and we do 
not discuss it here. Another important property is the 
superfluid density to be discussed next. 



E. Superfluid Density Tensor 

The superfluid density tensor is defined as a response 
function to phase twists of the order parameter—. In 
the approximation used in this paper, the temperature 
dependence of its components is given by 



i(k)s^(k) - r(k)a i e(k)a J c(k)] , (39) 



(38) 



where n(k) is the momentum distribution, Y(k) is the 
Yoshida function defined in the previous section, and di 
denotes the partial derivative with respect to k%. Gener- 
ally speaking, there are two components to the reduction 
of the superfluid density at non-zero temperatures, one 
coming from fermionic (quasi-particle excitations) and 
the other bosonic (collective modes) degrees of freedom. 
Here, we do not discuss the bosonic contribution, ex- 
cept to say that at low temperatures the dominant terms 
come from Goldstone modes associated with the phase 
of the order parameter—, which are undcrdamped in our 
case due to sub-critical Landau damping. Furthermore, 
in the case of tetragonal symmetry, Goldstone modes do 
not contribute to the off-diagonal component of the su- 
perfluid density, which is the main focus of the analysis 
that follows. 

We discuss first py for SHS p-wave superfluids. Plots 
of p l J — Pij/p™y X are shown in Fig.Elas a function of tem- 
perature for case (I), and three filling factors: iV c = 0.5 
(line), N c = 0.45 (hollow squares), and iV c = 0.4 (solid 
squares). The normalization constant is the maxi- 
mum value of the off-diagonal component in a square lat- 
tice. This maximum occurs at half-filling and T w 0.02T 
in a square lattice, where Tq is the critical temperature at 
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FIG. 6: Plots of reduced diagonal (p xx = p yv \ upper curves) 
and off-diagonal (p xy — pf. x ; lower curves) superfluid density 
tensor components pi? = Pij/p™y* versus temperature T/T c 
for SHS p-wave superfluids in tetragonal lattices (case (I)). 
N c — 0.5, 0.45 and 0.4 are shown with line, hollow and solid 
squares. 




FIG. 7: Plots of reduced diagonal (pf. x — p yy ; upper curves) 
and off-diagonal (p xy — p yx ; lower curves) superfluid density 
tensor components pi? — Pij/p™y* versus number of atoms per 
unit cell N c for SHS p-wave superfluids in tetragonal lattices 
(case (I)). Ti = 0.002T , T 2 = 0.02T and T 3 = 0.1T are 
shown with line, solid and hollow squares. Here To is the 
critical temperature at half-filling. 



half-filling. It is important to emphasize that square lat- 
tices [case (I)] have identical diagonal elements p xx = p yy 
due to tetragonal symmetry, but have nonzero p xy com- 
ponent as a result of the absence of reflection symmetry 
in the yz-plane (x — > —x) or the xz-plane (y — > —y) for 
the d- vector defined in Eq. iffijjl . In cases (II) and (III), 
reflection symmetry is restored and p xy vanishes identi- 
cally in the orthorhombic case for any temperature T. 

In the case of tetragonal symmetry [case (I)], p xy has 
a peak as a function of N c , but the diagonal components 
have a dip of the same size at exactly half-filling. In 
Fig. [7| we plot Pij/p^y* as a function of N c for three 
temperatures. 

Furthermore, notice that at T — the superfluid den- 



sity tensor is reduced to 

^ = 2^X>( k )W( k )> ( 4 °) 

which is just the integral of the momentum distribution 
weighted by the the curvature tensor of the dispersion 
£(k). Thus, p xy is zero at T = for any filling fac- 
tor. However, as T increases p xy increases and reaches 
approximately the same values as p xx and p yy . Further 
increase of T leads all tensor components of ptj to vanish 
at T = T c , as expected. 

In contrast, for s-wave superfluids, the off-diagonal 
tensor elements (p xy ,p yx ) are zero both in square (case 
(I)) and orthorhombic lattices (cases (II) and (III)). How- 
ever, similar to SHS p-wave superfluids, p xx = p yy in the 
square and p xx ^ p yy in the orthorhombic lattices. We 
do not plot p xx and p yy for the s-wave case, because our 
main interest here is the analysis of p xy which is strictly 
zero in this case. We note in passing that the main dif- 
ference between the diagonal components of s-wave and 
SHS p-wave superfluid densities is their temperature de- 
pendence. In the s-wave case at low temperatures, p xx 
and p yy deviate from their zero temperature value only 
by a small exponential correction, while in the SHS p- 
wave case the deviation is a power law. 

Measurements of the superfluid density in cold atoms 
might be performed through the rotation of the conden- 
sate in combination with experimental techniques that 
are sensitive to rotations. Through this rotational sen- 
sitivity it may be possible to extract the velocity fields 
and their correlations, as it is currently possible with 
positional (density) sensitive techniques^ 3 -. This experi- 
ment in cold atoms would be the analogue of the rotating 
bucket (or cylinder) experiments in liquid Helium 42 . 

In this section, we have discussed some of the exper- 
imentally relevant quantities of SHS p-wave superfluids 
and compared them with the s-wave superfluids in the 
weak coupling BCS regime. Furthermore, we presented 
signatures of the SHS p-wave superfluid state, which can 
be used to identify superfluidity in this limit. In the 
next section, we study some of these thermodynamic 
quantitites as a function of interaction strength from 
weak (BCS) to strong (BEC) coupling regimes at low fill- 
ing factors, and complete our analysis with a discussion 
of signatures of a possible BCS-to-BEC quantum phase 
transition for SHS p-wave superfluids. 

IV. SUPERFLUID PROPERTIES IN THE 
BCS-TO-BEC EVOLUTION 

At very low temperatures (T « 0), the saddle point 
order parameter and the number equations are approxi- 
matelly valid for all couplings from weak to strong cou- 
pling regimes, where a small correction to the number 
equation is negligible 38 . 

For low filling factors < N c < 0.5, the chemical po- 
tential fi = —\jx\ decreases as a function of interaction 
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strength Vo, x and crosses the bottom of the energy band 
(emin = —tx ~ ty) at some critical value of Vq, x . The 
decrease of fi is associated with the formation of bound 
particle pairs which are pulled out of the two-particle con- 
tinuum. Similarly, for high filling factors 0.5 < N c < 1, 
Ji = increases with increasing Vo, x and crosses the top 
of the energy band e max = t x + t y at some critical value of 
Vo,x- The increase of fi is associated with the formation 
of bound hole pairs which are pulled out of the two-hole 
continuum. However, exactly at half-filling N c — 0.5, jl is 
pinned at zero for any interaction strength due to perfect 
particle-hole symmetry. 

Notice that the situation in lattices is very different 
from the continuum, because of particle-hole symmetry. 
In the continuum the chemical potential is only pulled 
down below the bottom of the band since there is no 
energy upper bound for a parabolic dispersion. In the 
lattice discussed above, the chemical potential 

can move below (above) the bottom (top) of the band 
for filling factors below (above) half-filling. Notice in ad- 
dition, that for optical lattices of cold fermion atoms at 
high filling factors, superfluidity is associated with corre- 
lated motion of holes (atom voids), rather than particles 
(atoms). 



A. Phase Diagram 

Depending on the behavior of //, the T = phase di- 
agram can be plotted as a function of filling factor N c 
and interaction strength Vq >x . The BCS region, where 
\fl\ < t x + t y , corresponds to the weaker interaction re- 
gion for fixed density as shown in Fig. |SJ In the BCS 
region there is a well defined Fermi surface locus where 
£(k) = 0. The BEC region, where |/2| > t x + t y , corre- 
ponds to the stronger interaction region for fixed density 
as shown in Fig. |S| In the BEC region there is no Fermi 
surface locus, since £(k) can not be zero. 
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FIG. 8: Plots of critical interaction strength V r = Vo,x/Eo 
and V r — Vo )S /Eo versus number of atoms per unit cell N c 
for a) SHS p-wave, and b) s-wave superfluids, respectively at 
T = 0. Plots for a = 0,0.1 and 0.6 are shown with hollow 
squares, solid squares and line, respectively. 



In Fig.|H|wc also compare the phase diagrams of SHS p- 
wave and s-wave superfluids. In case (I) (hollow squares) 
of SHS p-wave superfluids, a wide region of filling factors 
around half-filling requires very strong critical Vo. x to 
reach the BEC regime. For a small anisotropy a = 0.1 
(solid squares) corresponding to case (II), this region nar- 
rows. However, further anisotropy a — 0.6 (line), this 
region widens again. For s-wave superfluids, in case (I) 
(hollow squares), the region around half- filling expands 
very little for small anisotropy a = 0.1 (solid squares) 
corresponding to case (II). Further anisotropy a = 0.6 
(line) expands the BCS region around half-filling, thus 
making it more difficult to reach the BEC regime at fixed 
filling factor. 

Notice in addition, that when the filling factor iV c — ► 
or N c — > 1 a finite interaction strength is necessary to 
evolve from the BCS to the BEC regime in both SHS 
p-wave and s-wave superfluids. This indicates that two- 
particle (or two-hole) bound states in a two-dimensional 
lattice require finite energy for the symmetries discussed. 
This is in contrast with the situation found in two dimen- 
sional continuum models, where a two-body bound state 
is found at arbitrarily small attractive interaction for 
s-wave. However, in two-dimensional continuum mod- 
els for the SHS p-wave case, the creation of a two-body 
bound state requires a finite interaction strength. 

Furthermore, there are two possible ways of investigat- 
ing the evolution from BCS to BEC regime in cold atom 
experiments. The first way is by changing the interac- 
tion strength for a fixed filling factor, while the second is 
by changing filling factor for a fixed interaction strength. 
Probably both ways are possible in cold atom experi- 
ments, where interaction strength and atom filling factor 
could be tuned independently 

In all physical properties to be discussed in the next 
sections, first we fix the filling factor to N c = 0.25 
(quarter-filling) and vary V r to cross the BCS-BEC 
boundary. But, in addition, we fix the interaction 
strength to V r = 6, and vary N c to cross the BCS-BEC 
boundary. We present results for tetragonal (case (I) 
with a = 0) and orthorhombic (case (II) with a = 0.1) 
lattices. Notice that cases (II) and (III) have similar be- 
havior, and differ only by a scale factor. Thus we do not 
present case (III) here. 



B. Quasi-particle Excitation Spectrum 

The quasi-particle excitation spectrum of SHS p-wave 
and s-wave superfluids during the BCS to BEC evolu- 
tion are very different because of the symmetry of the 
order parameter. As discussed in Sec. [nj the quasi- 
particle sp ectrum of SHS p -wave superfluids is given by 
E(k) = v/C 2 (k) + |A(k)| 2 with f (k) = e(k) - /L Since 
the cross-product ef*(k) x d(k) ^ is non-zero, it is 
expected that triplet superfluids have additional quasi- 
particle excitation branches. However, for the SHS p- 
wave state these branches do not enter the problem as 
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they are pushed to extremely high energies and are not 
accessible. For instance, when a THS p-wave superfluid is 
formed from a pseudo-spin 1 /2 system, there are two ac- 
cessible quasi-particle branches when time reversal sym- 
metry is broken. Thus, in the SHS p-wave there is only 
one quasi-particle energy branch. 

For the SHS p-wave case, notice that E(k) = when 
both £(k) = and A(k) = 0. This implies that 
the momentum space region of zero -E'(k) occurs when 
— t x cos(k x a x ) — t y cos(k y ay) — jl and Aq, x siri(k x a x ) + 
Ao,j, sm(k y a y ) — for chemical potentials inside the BCS 
region < t x +t y ). This means that in the BCS region 
the quasi-particle excitation spectrum is gapless. How- 
ever, E(k) never vanishes for chemical potentials inside 
the BEC region > t x + t y ), since £(k) can never be 
zero. However, the order parameter A(k) can still be 
zero. This implies that in the BEC region the quasi- 
particle excitation spectrum is fully gapped. Notice that 
this change in minimum excitation energy accompanies 
the existence or non-existence of the Fermi surface lo- 
cus £(k) = 0. Since quasi-particle excitations are evolv- 
ing from gapless BCS to gapped BEC regime as a func- 
tion of interaction strength for fixed filling factor, or as 
a function of filling factor for fixed interaction strength, 
the evolution between the BCS and BEC regimes is not 
smooth and a quantum phase transition occurs. As we 
shall see in Sec. IIVDI where ground state properties like 
the isothermal compressibility are calculated, this quan- 
tum phase transition is characterized by the chemical po- 
tential crossing either the bottom (jl c = — t. 
top {jl c = t x + t y ) of the energy band. 

In contrast, for the s-wave case, E(k) never vanishes 
and is always gapped in both the BCS and BEC regimes. 
This is a major difference between SHS p-wave (or more 
generally any nodal superfluids) and s-wave superfluids. 
In s-wave superfluids quasi-particle excitations are al- 
ways gapped and the evolution from the BCS to the 
BEC regime is smooth, so that no quantum phase tran- 
sition occurs. This gapless to gapped quantum phase 
transitions were first considered in the context of He 3 by 
Volovik 47 , and later in the context of high-T c supercon- 
ductivity^2iiS and atomic Fermi gasesii^i 3 -. 



C. Chemical Potential 

Plots of the reduced chemical potential p r = jl/Eo 
at low temperatures (T w 0) and quarter filling factor 
(N c = 0.25) are shown in Fig. as a function of reduced 
interaction V r — Vq :X /E and V r = Vo :S /E for SHS p- 
wave and s-wave superfluids, respectively. The tetrago- 
nal case (I) with a = and orthorhombic case (II) with 
a = 0.1 are shown with hollow and solid squares, respec- 
tively. In the case of SHS p-wave superfluids, p r is a 
continuous function of V r , but its slope with respect to 
V r has a cusp when jl crosses the bottom of the energy 
band; fx r = — 1 (a = 0: hollow squares) and /i r = —1.1 
(a = 0.1 :solid squares). However in the s-wave case, 



both fi r and the slope of p r with respect to V r are con- 
tinuous functions of V r . 
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FIG. 9: Plots of reduced chemical potential fi r = fi/Eo, 
versus reduced interaction strength Vr = Vo, x /Eo and Vr = 
Vo, 3 /Eo for a) SHS p-wave, and b) s-wave superfluids, respec- 
tively at T = and iV c = 0.25. Plots for a = and 0.1 are 
shown with hollow and solid squares, respectively. 

Plots of the reduced chemical potential p r = jl/Eo 
at low temperatures (T « 0) and constant interaction 
strength (V r = 6) are shown in Fig. ^5] as a function of 
filling factor iV c for both SHS p-wave and s-wave super- 
fluids. Here again the tetragonal case (I) with a = and 
orthorhombic case (II) with a = 0.1 are shown with hol- 
low and solid squares, respectively. In the case of SHS 
p-wave superfluids, p r is a continous function of N c , but 
its slope with respect to N c has a cusp when p, crosses 
the bottom of the energy band; fi r = — 1 (a = 0: hol- 
low squares) and \i r = —1-1 (a = 0.1: solid squares). In 
addition, p r changes curvature at p r = (N c = 0.5) at 
the same place, where the topology of the Fermi surface 
locus £(k) = changes from particle- like to hole- like. 
However, in the s-wave case, both (j, r and the slope of \i r 
with respect to N c are continuous functions of N c , and 
there is no change of slope when the topology of Fermi 
surface locus £(k) = changes. Thus, the slope change of 
/i r for the SHS p-wave case is directly related to a higher 
angular momentum pairing channel (£ ^ 0). 

This difference between SHS p-wave and s-wave super- 
fluids have great importance in describing the evolution 
from BCS to BEC regime. As discussed above, the major 
qualitative differences between the chemical potentials of 
SHS p-wave and s-wave cases are more clearly seen in 
the derivative djl/dN c , which is directly related to the 
atomic compressibility to be discussed next. 



D. Atomic Compressibility 

Plots of the reduced isothermal atomic compressibil- 
ity K r = kt(T)/kq at low temperatures (T w 0) (see 
Eg. I34f) and quarter filling factor (N c = 0.25) are shown 
in Fig. 1111 as a function of reduced interaction strength 
V r = Vq, x /Eq and V r = Vq, s /Eq for SHS p-wave and s- 
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FIG. 10: Plots of reduced chemical potential fi r = jl/Eo, 
versus filling factor N c for a) SHS p-wave, and b) s-wave su- 
perfluids at T = and V r = 6. Plots for a — and 0.1 are 
shown with hollow and solid squares, respectively. 



FIG. 12: Plots of NcK r = dN c /d^ r versus filling factor iV c 
for a) SHS p-wave, and b) s-wave superfluids at T = and 
K- = 6. Plots for a = and 0.1 are shown with hollow and 
solid squares, respectively. 



wave superfluids, respectively. The tetragonal case (I) 
with a = 0, and the orthorhombic case (II) with a = 0.1 
are shown with hollow and solid squares, respectively. 
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FIG. 11: Plots of reduced isothermal atomic compressibility 
Kr = kt/k-o versus reduced interaction strength V r = Vo, x /Eo 
and V r — Vo,s/Eq for a) SHS p-wave, and b) s-wave superflu- 
ids, respectively, at T — and N c = 0.25. Plots for a = 
and 0.1 are shown with hollow and solid squares, respectively. 

In the case of SHS p-wave superfluids, the evolution of 
K r from BCS to BEC regime is not smooth, and a cusp 
occurs at critical interaction strengths corresponding to 
fi r = —1 in case (I) (hollow squares) and [i r = —1.1 
in case (II) (solid squares). These cusps are associated 
with the appearance of a full gap in the quasi-particle 
excitation spectrum as the evolution from BCS to BEC 
takes place. In contrast, the evolution of smooth 
for s-wave superfluids, since the quasi-particle spectrum 
is always gapped during the BCS to BEC evolution. 

Plots of 2VgK r = dN c /d/i r at low temperatures (T w 0) 
and at fixed interaction strength (V r = 6) are shown in 
Fig. E| as a function of filling factor N c for both SHS 
p-wave and s-wave superfluids. The tetragonal case (I) 
with a = and orthorhombic case (II) with a — 0.1 are 
shown with hollow and solid squares, respectively. 

In the case of SHS p-wave superfluids, the evolution of 



N^K r — dN c /d[i r from the BCS to the BEC regime is 
not smooth. In both cases (I) (hollow squares) and (II) 
(solid squares) , cusps occur for corresponding filling fac- 
tors when /i r = ±1 and fx r = ±1.1, respectively. These 
cusps occur again as a result of the gapless to gapped 
transition in the quasi-particle excitation spectrum dis- 
cussed above. Notice that the cusps appearing at lower 
filling factors are associated with particle BCS-BEC tran- 
sitions, while the cusps appearing at higher filling factors 
are associated with hole BCS-BEC transitions. Further- 
more, N^n r has a sharp peak at half-filling in case (I), 
and two small peaks in case (II) which are symmetric 
around half-filling. These peaks occur inside the BCS re- 
gion of the phase diagram shown in Fig. [HI an d they arise 
due to the same reasons discussed in Sec. lIII Dl where the 
BCS regime is extensively analysed. In contrast, the evo- 
lution of N^K r as a function of N c is again smooth for 
s-wave superfluids during the BCS to BEC evolution. 

In cold Fermi gases the measurement of the isother- 
mal compressibility kt(T) at low temperatures is very 
hard, because most measurements are performed when 
traps are turned off and the gas expands. However, the 
gas expansion process is probably close to being isen- 
tropic and k$ is most likely measurable. As discussed 
in Sec. 



IIII Dl for T T c , the isentropic compressibility 
and the isothermal compressibility are essentially propor- 
tional, and their measurements can serve as an indicator 
of the quantum phase transition discussed in this section. 
Another important quantity that reflects such a transi- 
tion is the Cooper pair size to be discussed. 



E. Average Cooper Pair Size 

The average Cooper pair size in the saddle point ap- 
proximation is given by 

^ air = (^(k)|(-VDb(k))/(^(k)|^(k)) (41) 



14 



where </?(k) = A(k)/2£ , (k) plays the role of the Cooper 
pair wave function. This quantity reflects the average 
size of a Cooper pair and not the coherence length £ co h 
of the superfluid. Although these quantities are directly 
related in the BCS regime, they are very different in the 
BEC regime^. Here, however, we concentrate only on 
the analysis of £ pa ir- 

Plots of the reduced average Cooper pair size = 
^pair/a at low temperatures (T ss 0) and at quarter filling 
(N c = 0.25) are shown in Fig. ^31 The plots are shown as 
a function of reduced interaction strength V r — Vo, x /Eq 
and V r — Vo s /Eq for SHS p-wave and s-wave superfluids, 
respectively. The tetragonal case (I) with a = 0, and 
the orthorhombic case (II) with a = 0.1 are shown with 
hollow and solid squares, respectively. 
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FIG. 13: Plots of reduced average Cooper pair size £ r = 
Cpair/a versus reduced interaction strength V r = Vo, x /Eo and 
V r = Vo, s /Eo for a) SHS p-wave, and b) s-wave superfluids, 
respectively, at T = and iV c = 0.25. Plots for a = and 0.1 
are shown with hollow and solid squares, respectively. 



a. This is a manisfestation of higher angular momentum 
pairing and of the Pauli exclusion principle. 

In contrast, the evolution of £ r is smooth for s-wave 
superfluids, since the momentum distribution n(k) never 
vanishes at low k, and thus the distribution of pair sizes 
is well behaved. This monotonic decrease of £ pa ir is also 
a reflection of a quasi-particle spectrum that is always 
gapped during the BCS to BEC evolution and of an or- 
der parameter for zero angular momentum pairing. No- 
tice that, £ r is decreasing monotonically as a function of 
interaction V r for N c = 0.25. The limiting pair size for 
large V r and fixed N c is small in comparison to the lattice 
spacing a, since the Pauli principle does not forbid atoms 
of opposite pseudo-spins to be on the same lattice site. 

Next, we discuss the behavior of £ r as a function of fill- 
ing factor N c . Plots of £ r = £ pa i r /a at low temperatures 
(T s» 0) and at fixed interaction strength (Vo, x — 6J3o) 
are shown in Fig. [21 as a function of N c for both SHS 
p-wave and s-wave superfluids. The tetragonal case (I) 
with a — and orthorhombic case (II) with a = 0.1 are 
again shown with hollow and solid squares, respectively. 





In the case of SHS p-wave superfluids, the evolution of 
£ r from BCS to BEC regime is not smooth, and singular 
behaviors occur at critical interaction strengths corre- 
sponding to fj, r = —1 in case (I) (hollow squares) and 
/i r = —1.1 in case (II) (solid squares). These singular be- 
haviors are associated with the appearance of a full gap 
in the quasi-particle excitation spectrum and with pair- 
ing in a non-zero angular momentum as the evolution 
from BCS to BEC takes place. 

This singular behavior can be understood in terms 
of the vanishing of the momentum distribution when 
jl smaller (larger) than the bottom (top) of the band 
for fixed filling factor N c . For instance, in the case of 
N c = 0.25 shown in Fig. ^5] , when p, falls slightly below 
the bottom of the band, suddenly u(k) vanishes in the 
neighborhood of k = which leads to unbound pairs of 
atoms with (ki = k = 0; f) and (k.2 = — k = 0; f), and 
thus a rapid increase in the average pair size £ pa ; r for 
interaction strengths below or above the critical values. 
Beyond the critical interaction strength, £ pa i r decreases 
monotonically for increasing V r and converges asymptot- 
ically (when V r — ► oo and N c < 0.5) to a finite value 
which is larger but of the order of the lattice spacing 



FIG. 14: Plots of reduced average Cooper pair size £ r = 
£pair/a versus filling factor iV c for a) SHS p-wave, and b) s- 
wave superfluids at T — and V r = 6. Plots for a — and 
0.1 are shown with hollow and solid squares, respectively. 

In the case of SHS p-wave superfluids, the evolution of 
£ r from BCS to BEC regime is not smooth. In both cases 
(I) (hollow squares) and (II) (solid squares), singulari- 
ties occur for corresponding filling factors when ji r = ±1 
and fi r = ±1.1, respectively. These singularities occur 
again as a result of the gapless to gapped transition in 
the quasi-particle excitation spectrum discussed above, 
as well as the symmetry of the order parameter, which 
strongly modifies the nature of the pair wave function 
<ys(k) = A(k)/22?(k). Notice that £ r has a singularity at 
half- filling (BCS region) in case (I), and two other sin- 
gularities symmetric around half-filling (BCS region) in 
case (II) . As the interaction is further increased, the low 
and high filling singularities (BCS-BEC boundaries) shift 
towards half-filling and merge with the half-filling singu- 
larity in case (I) (see also phase diagram in Fig.[SJ). This 
situation is analogous to the ci-wave lattice cas o 50 ' 51 . As 
the interaction is further increased, similar behavior oc- 
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curs in case (II). Thus, at large interaction strengths the 
BCS (BEC) region shrinks (expands), and the character- 
istic value of £ pair is of the order of the lattice spacing 
a except when iV c is very close to half-hlling, where the 
BCS-BEC phase boundary is located. 

In contrast, the evolution of £ r as a function of N c is 
again smooth for s-wave supcrfluids during the BCS to 
BEC evolution (see Fig. 114b). For instance, below half- 
hlling £ r decreases as a function of increasing N c (increas- 
ing DOS), while the pair wave function <p(k) continues 
having the same qualitative behavior at low momenta. 
Above half- filling, analogous discussions hold for holes. 
Furthermore, £ pa ir decreases monotonically with increas- 
ing interaction for any filling factor N c , indicating that 
for very large interaction strengths £ pa i r — > for any N c . 

Thus far, we have investigated only momentum aver- 
aged quantities, but next we discuss the momentum dis- 
tribution, which could be measured in cold Fermi gases. 



Momentum Distribution 



The zero temperature momentum distribution 



n(k) = 



1 



1 - 



(42) 



is discussed in this section in the weak coupling BCS and 
strong coupling BEC regimes for both SHS p-wave and 
s-wave superfluids in a lattice. We discuss here only the 
case of particle superfluidity with N c — 0.25. In the BCS 
regime, we choose V r corresponding to (x r = —0.5 for 
case (I) (a = 0) and /z r = —0.55 for case (II) (a = 0.1). 
Similarly in the BEC regime, we choose V r corresponding 
to n r = — 2 for case (I) and fi r = —2.2 for case (II). 

In the case of SHS p-wave superfluids, the momentum 
distribution has a major rearrangement in k-space with 
increasing interaction strength as can be seen in Fig. IT51 
This rearrangement is very dramatic when // crosses the 
bottom of the energy band /2 = fl c . This is a direct 
consequence of the change of the quasi-particle excitation 
spectrum from gapless in the BCS regime < (2 + a)t) 
to fully gapped in the BEC regime (|/2| > (2 + a)t). 

For case (I), plots of n(k) (in the first Brillouin zone) 
for BCS and BEC regimes are shown in Fig. 115b and 
Fig. 115b. respectively. Here the interaction strengths cor- 
respond to fj, r = —0.5 for the BCS regime, and \i r = —2 
for the BEC regime. The boundary line between the two 
regimes occurs when fi r = —1. Notice that n(k) is sym- 
metric around the lines k y a y = k 

x&x and kyCLy — k x cc x 
in both regimes as a reflection of the symmetry properties 
of £(k) and |A(k)|. 

While the line k y a y = —k x a x has the highest mo- 
mentum distribution in the BCS regime, n(k) vanishes 
along k y a y = —k x a x in BEC regime. As the interac- 
tion strength increases, two-particle states with opposite 
momenta are taken out of the two-particle continuum 
into two-particle bound states with zero center of mass 
momentum. As more of these tightly bound states are 
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FIG. 15: (Color online) Plots of momentum distribution n(k) 
versus k x a x and k y a y for an SHS p-wave superfluid in case (I) 
in a) BCS (fi r = —0.5) and b) BEC (fi r = —2) regimes. In 
addition plots for case (II) in c) BCS (p r = —0.55) and d) 
BEC (fi r = —2.2) regimes at T = 0. Here the brighter the 
region, the higher the value of momentum distribution. 



formed the large momentum distribution in the vicinity 
of k = splits into two peaks around finite momentum 
values reflecting the t = 1 value of the angular momen- 
tum associated with these p-wave tightly bound states 
in the BEC regime. For instance, in the case of d-wave 
(£ = 2) superfluids the momentum distribution in the 
BCS regime is centered around k = and splits into 
four peaks around finite momentum values reflecting the 
£ = 2 value of the angular momentum associated with 
these tightly bound states in the BEC regime^ 

For case (II), plots of ?i(k) for BCS and BEC regimes 
are shown in Fig. 115b and Fig. I15H . respectively. The 
figures shown correspond to interaction strengths asso- 
ciated to /!,. = —0.55, and fi r = —2.2 in the BCS and 
BEC regimes, respectively. The boundary line between 
the two regimes occurs when /i r = —1.1. Notice that, 
rt(k) is symmetric around k x a x = and k y a y = in 
both regimes. While the highest values of rt(k) occur 
near k = in the BCS regime, n(k) vanishes near k = 
and develop two maxima for finite value of k x in the BEC 
regime. 

In contrast, the evolution of the momentum distribu- 
tion from BCS to BEC regime is smooth for s-wave su- 
perfluids as shown in Fig. ^| for case (I). In case (II) 
the results are very similar to case (I), except for the 
expected anisotropy and we do not present them here. 
Notice that for s-wave (£ = 0) superfluids, the momen- 
tum distribution in the BCS regime is centered around 
k = 0, and remains centered around k = even in the 
BEC regime. This reflects the £ = value of the an- 



16 



(a) (b) 




-3-2-10123 -3-2-10123 



k x a x k x a x 

FIG. 16: (Color online) Plots of momentum distribution n(k) 
versus k x a x and k y a y for an s-wave superfluid in case (I) 
in a) BCS (fx r = —0.5) and b) BEC (/i r = —2) regimes at 
T — 0. Here the brighter the region, the higher the value of 
meomentum distribution. 

gular momentum associated with s-wave tightly bound 
states in the BEC regime. Here, the essential qualitative 
difference in the BCS to BEC regimes is that as the inter- 
action gets stronger n(k) broadens due to the formation 
of tightly bound states. 



V. CONCLUSIONS 

In summary, we considered SHS p-wave pairing of sin- 
gle hyperfine state and THS s-wave pairing of two hyper- 
fine state Fermi gases in quasi-two-dimensional optical 
lattices. We focused mainly on superfluid SHS p-wave 
(triplet) states that break time-reversal, spin and orbital 
symmetries, but preserve total spin-orbit symmetry. 

Since the paper was divided into two parts, we present 
our conclusions for each of these parts separately. In the 
first part, we analysed superfluid properties of SHS p- 
wave and THS s-wave symmetries in the strictly weak 
coupling BCS regime. There, we calculate the order pa- 
rameter, chemical potential, critical temperature, atomic 
compressibility, and superfluid density as a function of 
filling factor for tetragonal and orthorhombic optical lat- 
tices. 

We found that, for SHS p-wave superfluids, the crit- 
ical temperatures in tetragonal and orthorhombic opti- 
cal lattices are considerably higher than in continuum 
model predictions, and therefore, experimentally attain- 
able. In particular, we showed that a small anisotropy 
in the transfer energies increases the effective density of 
fcrmionic states of an SHS p-wave superfluid considerably 
around half-filling which leads to higher critical tempera- 
tures. In contrast, a small anisotropy decreases the den- 
sity of states around half-filling, and therefore, the criti- 
cal temperatures for an THS s-wave superfluid. We also 
noticed that further anisotropy leads to a larger decrease 
in critical temperatures for THS s-wave than for SHS p- 
wave superfluids. Therefore, we concluded that a small 
anisotropy favours SHS p-wave pairing near half-filling, 
and that the critical temperatures for SHS p-wave super- 
fluids are comparable and even higher than THS s-wave 



case for similar interaction parameters. 

For SHS p-wave superfluids, we found a peak in the 
isothermal atomic compressibility at low temperatures, 
exactly at half-filling of tetragonal lattices. This peak 
splits into two smaller peaks in the orthorhombic case. 
These peaks reflect the SHS p-w&ve structure of the order 
parameter at low temperatures and they are related to 
the nodes (zeros) of the quasi-particle energy spectrum. 
We also showed that the atomic compressibility peaks 
decrease in size as the critical temperature is approached 
from below. The peaks turn into humps at T = T c . In 
contrast, for THS s-wave superfluids, we showed that the 
atomic compressibility does not show a peak structure at 
low temperatures since the quasi-particle energy spec- 
trum is always gapped, and that the compressibility is 
largely temperature independent. 

We also discussed the superfluid density tensor for SHS 
p-wave and THS s-wave systems. In the SHS p-wave 
case, we concluded that in orthorhombic lattices, the off- 
diagonal component p xy of the superfluid density tensor 
vanishes identically, while the diagonal components p xx 
and p vv are different. However, in tetragonal lattices, 
we showed that p xy ^ 0, while p xx = p yy . In contrast, 
for THS s-wave superfluids, p xy — and p xx ^ p yy in 
the orthorhombic, while p xy = and p xx — p yy in the 
tetragonal lattices. Therefore, the presence of non-zero 
p xy in the square lattices is a key signature of our exotic 
SHS p-wave triplet state. 

In the second part of the manuscript, we analysed su- 
perfluid properties of SHS p-wave and THS s-wave su- 
perfluids in the evolution from BCS to BEC regime at 
low temperatures (T « 0). We discussed the changes in 
the quasi-particle excitation spectrum, chemical poten- 
tial, atomic compressibility, momentum distribution and 
Cooper pair size as a function of filling factor and inter- 
action strength for tetragonal and orthorhombic optical 
lattices. We found major differences between SHS p-wave 
and THS s-wave superfluids in the evolution from BCS 
to BEC regime. 

In the case of SHS p-wave superfluids, the quasi- 
particle excitation spectrum, chemical potential, atomic 
compressibility, Cooper pair size and momentum distri- 
bution are not smooth functions of filling factor or in- 
teraction strength. In particular, the singular behavior 
of the atomic compressibility suggests the existence of a 
quantum phase transition when the chemical potential 
crosses the bottom or the top of the energy band. This 
transition is associated with the change in quasi-particle 
excitation spectrum and the higher angular momentum 
nature (t = 1) of SHS p-wave superfluids, which are gap- 
less in the BCS regime, but fully gapped in the BEC 
regime. 

In contrast, for THS s-wave superfluids, the quasi- 
particle excitation spectrum, chemical potential, atomic 
compressibility, Cooper pair size and momentum distri- 
bution are smooth functions of filling factor or interac- 
tion strength. In particular, the smooth behavior of the 
atomic compressibility suggests only a crossover when 
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the chemical potential crosses the bottom or top of the 
energy band. This crossover is associated with a fully 
gapped quasi-particle excitation spectrum and pairing at 
the zero angular momentum channel in both BCS and 
BEC regimes. 

These major differences between THS s-wave and SHS 
p-wave superfiuids in an optical lattice suggest that SHS 
p-wave cold atoms are much richer than THS s-wave cold 
atoms, and thus these systems may provide a new exper- 



imental direction in this field. 
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